Electron spin decoherence in diluted magnetic quantum wells 
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We study electron spin dynamics in diluted magnetic quantum wells. The electrons are coupled by 
exchange interaction with randomly distributed magnetic ions polarized by magnetic field B. This 
coupling leads to both spin relaxation and spin decoherence. We demonstrate that even very small 
spatial fluctuations of quantum well width dramatically increase rate of decoherence. Depending 
on the strength of exchange interaction and the amplitude of the fluctuations the decoherence can 
be homogeneous or inhomogeneous. In the homogeneous regime, the transverse (with respect to B) 
component of electron spin decays on the short time scale exponentialy while the long-time spin 
dynamics is non-exponential demonstrating long-lived power law tail. In the inhomogeneous case, 
the transverse spin component decays exponentially with the exponent quadratic in time. 
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Effective manipulation of spin degree of free- 
dom in a semiconductor device by external elec- 
tric and magnetic fields is one of the primary 
goals of spintronicsjA A possible way to increase 
coupling of electron spin to external magnetic 
field is to use semiconductor-based materials, 
which incorporate magnetic elements. Such ma- 
terials are called diluted magnetic semiconduc- 
tors (DMS). The most common DMS are II- 
VI and IV- VI compounds with magnetic impu- 
rities (usually, Mn), and also III-V crystals^ 
like Gai_ x Mn x As (typically, x < 10%). The 
DMS combine magnetic and semiconductor prop- 
erties in a single material. This offers a large 
prospect for applications. In particular, DMS 
are considered to be the most promising can- 
didates in creating room-temperature ferromag- 
netic systems, which can be easily manipulated 
as semiconductors. Semiconductor heterostruc- 
tures doped by magnetic impurities have al- 
ready demonstrated exciting physical phenom- 
ena specific for magnetic systems: coherent spin 
excitations*^ magnetic polaron formation^ fer- 
romagnetic hole alignment^ etc. 

Many remarkable features of DMS, such as the 
large Zeeman splitting of the electronic bands 
and the giant Faraday rotation, are induced by 
the exchange interaction between the localized 
electrons on d-shells of Mn ions and delocalized 
band carriers (see Ref . [H for review) . This inter- 
action is also responsible for the collective nature 
of spin excitations in DMS, arising of novel collec- 
tive modesj2iiLIi, anticrossing of the electron and 
ion spin precession frequencies spin suscepti- 
bility enhancement ^ etc. 

The fluctuations of exchange field around the 



average value lead to electron spin relaxation 
and spin decoherence with characteristic times 
T\ and T2, respectively. Recently, the decoher- 
ence rate was measured in Cdi_ x Mn x Te two- 
dimensional (2D) structure.-^ The observed T2 
was found to be at least one order magnitude 
shorter than the decoherence time predicted the- 
oretically in Ref. [HI, where fluctuations of ex- 
change field were linked to delta-correlated fluc- 
tuations of magnetic ion concentration. In this 
paper, we demonstrate that even very small fluc- 
tuations of quantum well width dramatically in- 
crease rate of decoherence. Our estimates show 
that this mechanism can explain short decoher- 
ence time observed in Ref. [RJ 

We consider the 2D degenerate electron gas 
located in the {x, y) plane interacting with the 
magnetic ions randomly distributed with average 
2D concentration rij, and doping profile f(z) in 
the growth direction (J f(z)dz — 1). The system 
is placed into the magnetic field which we assume 
to be parallel to the well plane (B || e^) as it was 
the case in Ref. The field leads to Zeeman 
splitting of both electron and ion spin levels with 
energies huj e and hwj, respectively. The Hamil- 
tonian of the system is given by 



H = H e + Hj + H, 



Je • 



H e = h U(r) + hhj e s x , 

2m 

Hj = hojj^j kx , 

k 

H Je = asJ2^kS(r - R k )\^(z k )\ 
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Here H P is the Hamiltonian of an electron in 



the external magnetic field, r = (x, y) and p = 
iPxiPy) are respectively electron in-plane posi- 
tion vector and momentum, U(r) is the random 
impurity potential, which we assume to be short- 
range, Hj is the Hamiltonian of the ions and 
Hj e represents the exchange interaction between 
spin s of an electron placed at the lowest level in 
the well (with the wave function 9(z)) and the 
spins 3k of the ions located at points (Rk,Zk)- 
The strength of the interaction is characterized 
by constant a. 

It is convinient to rewrite Hj e as 



H, 



(H Je ) + SH 



Je, 



(5) 



where the angular brackets mean averaging over 
ion positions and thermal averaging and SHj e de- 
scribes fluctuations of exchange field. The term 
(Hj e ) leads to significant renormalization of the 
electron spin precession frequency 



*(Jx)nj 



f(z)\^(z)\ 2 dz, (6) 



Typically, u> ui e , 11,13 which means that the 
electron spin precession frequency is mostly de- 
termined by the effective magnetic field, created 
by polarized ions. 

The fluctuations of the exchage field arising 
due to random distribution of (Rk, Zk) lead both 
to electron spin relaxation and to decoherence. 
The analysis of these processes shows^ that in 
2D case the longitudinal and transverse (with re- 
spect to B) components of the electron spin de- 
cay exponentially with the characteristic times 
T[ ~ T' 2 ~ 2a 2 H 3 /3a 2 mnj. Comparing the re- 
sults of Ref. [li] with the recent experimental 
data^ one can see that experimentally observed 
decoherence time Ti is much shorter (by about an 
order of magnitude) than which implies that 
delta-correlated density fluctuations can not pro- 
vide sufficient fluctuations of the exchange field. 

Below we demonstrate that spatial fluctuations 
of quantum well width can (at certain condi- 
tions) lead to significantly shorter T2. Physi- 
cally, this happens because such fluctuations in- 
duce the long-range fluctuations of the effective 
magnetic field acting on the electron spin. 

First, we notice that the effective magnetic 
field induced by exchange interaction depends on 
the doping profile f(z) and quantum well ge- 
ometry (see Eq. ([6])). Expressions for T[ and 
T 2 were derived^ for the case of homogeneous 
distribution of magnetic ion and infinitely deep 
quantum well, having constant width ag. Let us 



now assume, for example, that the quantum well 
is infinitely deep and magnetic ions concentrate 
close to the center of the well (these assump- 
tions more or less correspond to the experimen- 
tal situation-^). Assuming also that well width 
slightly fluctuates, a(r) = ao + Sa(r) we find that 
spin precession frequency becomes r— dependent: 

/ \ 2ct(J x )nj . . . . 

w^w(r)=w e H = uj a + 6u(t), (7) 

na(r) 

where ujq = u> e + 2a(J x )nj/Hao, 5u)(r) ~ 
—2a(J x )nj5a(r) /ha^. We assume the fluctua- 
tions to be Gaussian with a spatial scale d: 



{5cj(r)6uj(r')) = 2u 2 X 



r — r 



(8) 



Here u 2 = 2a 2 {J x ) 2 n 2 (Sa 2 (r))/h 2 aQ is the am- 
plitude of the flucuations and x(£) i s a dimen- 
sionless function (x(0) = 1, %(£) — > 0, for £ s> 
1). As seen from Eq. ©, (5uj(r)6u>(r')) - (J x ) 2 , 
so that at low magnetic fields the fluctuations 
vanish and decoherence is due to the mecha- 
nism proposed in Ref. Il4l . However, at rela- 
tively strong fields corresponding to almost full 
polarization of ions (this was the case in the 
experiment^) the fluctuations are sufficiently 
large and dominate the decoherence. 

We will consider time evolution of the spin ex- 
citations concentrated near Fermi surface. The 
total transverse spin (per unit area) can be writ- 



ten as S+(t) 



iS z = e lUot S(t), where 



S(t) = J s + (r,ip,t)dipd 2 r/fl is slowly decaying 
amplitude, Q is the sample area and s + = s y +is z 
is the spin density that obeys the following qua- 
siclassical kinetic equatio n 15 ' 16 

2l± +VFn ^±- i s UJ (r)s + = St{s + }, (9) 

where tp is the velocity angle, vf is the Fermi 
velocity n = e^, cos (p+e y sin ip, and St{s + } is the 
collision integral describing the elastic scattering 
on the impurity potential U(r) with the mean 
free path I. The quasiclassical approach based 
on Eq. © is valid provided that d 3> h/mvp and 
hhj(r) <C rim p. Below we will solve this equation 
with the initial condition s+(r, ip, 0) = I. 

The spin decoherence can be homogeneous or 
inhomogeneous depending on the parameter w*r, 
where r is a characteristic time required for elec- 
tron to travel the distance of the order of d. For 
d <^ I, this time is given by r = d/vp, while for 
d > /, r = d 2 /AD ~ d 2 /lv F , where D = vpl/2 is 
the diffusion coefficient. 



For ujifT 3> 1, electron spins in different cor- 
relation regions rotate independently with local 
frequencies. Hence, the decoherence is inhomo- 
geneous and the transverse spin decays as 



S(t) 



(10) 



In the opposite case, uj*t <C 1, electron is visit- 
ing many correlated regions during decoherence 
time and the decoherence is homogeneous. First, 
we assume d <C I. If the inequality T^vf <C I 
is also satisfied, the electron motion on the time 
scale on the order of is ballistic. Neglecting 
the collision integral in Eq. ©, we obtain 

S(t) = ^/oM^nr)^ ^ 
— e ~ u l fo fa x(vF\T 1 -T 2 \/d)dT 1 dT 2 _ 

For t ^> d/vp, Eq. (fTT|) becomes 



S(t) 



- e ~t/T 2 



1 2wtd 



T2 VF 



X(0df- (12) 



where 



K-T't-v = 



r — r 



G D (r-r',t-t'), 



and Gd = (d/dt — DA) -1 is the Green function 
of the diffusion equation. Solving Eq. (|14l) with 
the initial condition s+(r, ip, 0) = 1, we find 



S(t) 



2tt 



- 4W + 2u; 2 d 2 f 00 7 ^L^i' 

* JO Dq z —iuj 2tt 



(15) 

where x is the Fourier transform of \- The 



It turns out that Eq. p2|) is also valid for T 2 vp ^> 
I. To see this, one can iterate Eq. ([9]) with 
respect to small 8u> and decouple correlations, 
which yields 1/T 2 = J dr' J rff(fo(r)G(r-r',i- 
i')&j(r')). Here G(r,t) is the Green function of 
Eq. © with d(jj(r) = (Green function of the 
Boltzmann equation) integrated over initial ve- 
locity directions and averaged over final velocity 
directions. This function can be presented as a 
sum over processes with different number of col- 
lisions: G(r,t) — (5(r — VFnt)) v exp(—tvF/l) + 
SiV=i CA r ( r i *)) where the first term is the ballis- 
tic contribution, which gives Eq. (|T2"j) . Relative 
contribution of other terms to I/T2 (compared 
to Eq. (|T2"jO may be shown to be on the order of 
parameter (d/l) ln(l/a;*T) <C 1^ 

For rf s> i and w»r = u}*d 2 /4D <C 1, the de- 
coherence is also homogeneous. In this case the 
electron motion is diffusive and one can by stan- 
dard means reduce Eq. (|9]) to diffusion equation: 



ds + 
~dt 



DAs + - iScu(r)s + = 0. 



(13) 



Iterating Eq. (p~3|) with respect to 6u> and decou- 
pling correlations, we obtain 



ds + 
~~dt 



DAt 



, t _ t ,s+(r',O = 0, 
(14) 



FIG. 1: Decay of the transverse spin amplitude in 
the diffusive case 

integrand in Eq. (|15|) has two poles and 
a branch cut along the negative imaginary 
axis. The poles are at the points W1.2 ~ 
(2x(0)uIt/tt) [±7r - 2£ln(l/Wr)] , where t = 
d 2 /4^,x(0)-2^/ oo ^(^. 




FIG. 2: Decoherence rate as a function of the correla- 
tion radius of the fluctuations. For d > do mechanism 
related to fluctuations of the well width dominates 

The imaginary parts of the poles are large com- 
pared to the real ones, so that the poles contri- 
bution is approximately given by 



The result looks 18 



S{t) w exp 



4x(0)w2rt 



hi 



(16) 



This contribution dominates at short times, 
w 2 ri <C 1. For uj 2 rt 3> 1, the main contribution 



is due to branch cut, yielding 
S(t) « - 



1 



2x(0) wjrt In 2 (t/r) ' 



(17) 



The contribution of the branch cut is negative. 
From Eqs. (|16l) and (fT7|) we see that the ampli- 
tude S(t) changes sign as shown in Fig. 1. We 
also see that a long-lived power-law tail appears 
in the transverse spin polarization. 

Above we assumed that B is parallel to the 
well plane and neglected the effect of the held on 
the orbital motion. The results are also valid for 
B || e z provided that R c » I (R c is the cyclotron 
radius). Our calculations can be easily general- 
ized for the opposite case, R c -C I. In particu- 
lar, in the homogeneous ballistic regime, under 
the assumptions w*r <C 1, d <C R c <C l,R c 
T2VF <C I, the transverse spin can be calculated 
in analogy with Eq. (|11[) by averaging of decoher- 
ence action calculated along ballistic trajectory. 



S(t) 



T 2 



ttR c 



x(m. (is) 



From equations derived above one can see that 
increasing correlation radius d decreases T2 both 
for d <C I and for d ^> I. The maximal value of 
I/T2 is on the order of w*. One can show that 
ui^T^ 3> 1 for typical values of parameters^ 
which implies that suggested mechanism domi- 
nates already at small d (d > do, see Fig. 2) 
and might be responsible for short values of de- 
coherence time observed in the experiment^ 
The parameter do is estimated as follows: do ~ 
hpa? /njSciQ, for I ^> d and do ~ ^/kpa^ljnjSa^, 
for I <C d regime. Taking Sag — yj (5a 2 (r)) on 
the order of the lattice constant, for typical val- 
ues of experimental parameters^ we hnd that 
do ~ lOnm. 
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